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1 Introduction 

In [3] and [6] (see also [4] ) , we present some definite integral and infinite series addition theorems 
which arise from expanding fundamental solutions of elliptic equations on M'^ in axisymmetric 
coordinate systems which separate Laplace's equation. We utilize orthogonality and integral 
transforms to obtain new definite integrals from some of these addition theorems. 

2 Definite integrals from integral transforms 

2.1 Application of Hankel's transform 

We use the following result where for x G (0, oo) we define 
F{r±0) := lim F{x); 

see [15, p. 456]: 

Theorem 1. Let F : (0, oo) C be such that 



roo 

/ \/lr\F{x)\dx < OO, 
Jo 



and let u > — Then 



/"OO /"OO 

-{F{r + Q)+F{r -{)))= j uJ^{ur) xF{x).J^{ux) dx du (2) 
2 Jo Jo 

provided that the positive number r lies inside an interval in which F{x) has finite variation. 

*This paper is a contribution to the Special Issue "Superintegrability, Exact Solvability, and Special Functions" 
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As an illustration for the method of integral transforms, we give the following example. 
According to [15, (13.22.2)] (see also [7, (6.612.3)]), we have for Re a > 0, 6, c> 0, Re i/ > 
then 



oo 



where Ju : C \ (— oo,0] — ?■ C, for order v £ C is the Bessel function of the first kind defined 
in [12, (10.2.2)] and : C \ (— oo, 1] — )• C for + ^ — N with degree v and order ^, is the 
associated Legendre function of the second kind defined in [12, (14.3.7), § 14.21]. The Legendre 
function of the second kind Qi, : C\ (— oo, 1] — )• C for ^ — N is defined in terms of the zero-order 
associated Legendre function of the second kind Quiz) := Q^i^)- If we apply Theorem 1 to the 
function F : (0, oo) — )• C defined by 

F{k) := "^e-'^Mkc), (4) 

then condition (1) is satisfied. If we use (4) in (2) then we obtain the following result. If 
Rea>0, c>0, Rei^> — |, then 

.2 I f,2 I „2 



which is actually given in [13, (2.18.8.11)]. 

Hardy [8, (33.16)] derives an interesting extension of (3) (see also [15, p. 389] and [2, p. 17]). 
We apply the Whipple formula [12, (14.9.17), § 14.21] to Hardy's extension to obtain 



ke'''''Jy{kb)Ju{kc)dk 

= Z2« 01 ( a' + b' + c^ \ 

^ V6c (a2 + (6 + c)^f' (a2 + (6 - c^f' ^"-^'^ \ 2bc J ' ^ > 

for Re a > 0, 6, c > 0, Re > —1. It is mentioned in [15, p. 389] that (5) can be derived from (3) 
by differentiation with respect to a. Using this integral and Theorem 1, we prove the following 
theorem. 

Theorem 2. Let Re a > 0, c > 0, Re i/ > -1. Then 

L {a^ + {b + cYf'\a- + {b-cn'^^^-'/A 26c j^^^^" 2a ^ ^^^^^^^ 
Proof. By applying Theorem 1 to the function F : (0, oo) — t- C defined by 

F{k) ■■= -"^e-'^Mkc), 
la 

using (5), we obtain the desired result. ■ 

Now, we give another example of how an integral expansion for a fundamental solution of 
Laplace's equation on in parabolic coordinates can be used to prove a new definite integral. 

Theorem 3. Let m G Nq, A' G (0,oo), /i, ^' G (0,oo), fi ^ fi' , k £ (0,oo). Then 

/•oo 

I Qm—1 

/2ix) Jm{kX)VXdX = 2n^/X^Jl4J/Jrnik\')LrnikfJ,<:)Krnik|I>), 
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where 

^ ~ 8XX'^^^' ^ 

r min r /-i 

and := 

Proof. We apply Theorem 1 to the function F : (0, oo) — )• C defined by 

F{k) := 2TTy^X' nfi' JmikX')ImikfJ-<)Km{k^i>), 

where 1^ : C \ (— oo, 0] — )• C, for z/ € C, is the modified Bessel function of the first kind defined 
in [12, (10.25.2)] and Ki, : C\ (— oo,0] — ^ C, for G C, is the modified Bessel function of the 
second kind defined in [12, (10.27.4)]. Again, we see that condition (1) is satisfied. We obtain 
the desired result from [6, (41)], namely, for A S (0,oo), 

Jjn{kX)Jm{kX')Im{kn<)Km{k^i>)kdk= ^"^ ■ 

27rvAA'/i/i' 

2.2 Application of Fourier cosine transform 

Theorem 4. Let a,b ^ (0, oo) with b < a, A; G (0,oo), Rei^ > Then 

Qv-l/2 ( - 

'0 

Proof. According to [7, (6.672.4)] we have the integral relation 

1 ^ fa^ + b'^ + z^^ 



f 

Jo 



/ a2 _^ f,2 _^ ^2 \ 

Qv-\j2 ( j cos(/c2:) dz = TTV ab Iy{kh)Ky{ka) . (6) 



1 / 

Iy{kh)Ky{ka) cos{kz)dk = — -=Q^_ij2 I 

2v ah \ 



2ab 

where a,b & (0, oo) with b < a, z > 0, Re u > We obtain the desired result from Theorem 1 
with F : (0, oo) — )■ C defined such that 

F{k) := T^\I^Iu{kh)K^{ka) 

V nj 

and = — |. Furthermore, if one makes the replacement z i— t- z/{^/2a), k i— t- \f2ka in [7, 
(7.162.6)], namely 



^ <3i/-i/2(l + z^) cos{kz)dz = -^Ijy (^-^ Ky ^-^^ 



where k G (0,oo) and Rei/ > — ^, then we see that (6) holds also for any a,b £ (0,oo) with 
a = b. m 

3 Definite integrals from orthogonality relations 

3.1 Degree orthogonality for associated Legendre functions 
with integer degree and order 

We take advantage of the degree orthogonality relation for the Ferrers function of the first kind 
with integer degree and order, namely (cf. [7, (7.112.1)]) 

/ F^icose)P^,{cose)sm9d9= ) ^ j ^^,^., (7) 
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where m,n,n' G Nq, and m < n, m < n' . We are using the associated Legendre function of the 
first kind (on-the-cut), '■ (— li 1) — ^ C, for u, fj, £ C, the Ferrers function of the first kind, 
which is defined in [12, (14.3.1)]. 

The following estimates for the Ferrers function of the first kind with integer degree and 
order will be useful. If 6* G [0, tt] and m, n e No then [14, § 5.3, (19)] 

|P-(cose)|<^^^^ (8) 
n! 

and if (9 G (0,7r) then [10, p. 203] 

|p-(cos^)| < 2 ^^+,""^' (vrn)-!/^ (,s,0)m+i/2^ 
n! 

If /i E C, ^ > are fixed and < z/ — )• +oo, we also have the following asymptotic formulas for 
the associated Legendre functions 

P,^(coshe) = (2^sinh0-^/^ ^^;,^^|.^^ e(-+^)^(l + O(^-^)), (10) 
P,^(isinhO = (27rcoshe)-V2 ^('' + /^ + ^) e('-+l)^+W2(i^o(^-i))^ (12) 

r(^' + 2) 

These asymptotic formulae follow from representations of Legendre functions by Gauss hyper- 
geometric functions; see [1, (8.1.1), (8.10.4-5)]. 

Theorem 5. Let n,m £ No, with n > m, u £ C \ {2m, 2m + 2, 2m + 4, . . r, r' E (0, 00), 
r / r', e' G (0,7r). Then 

[ ix' - i)(''+^)/^Qi\72^/^x)p;r(cos^)(sin^)('^+^)/2d0 

where 

X = ^ , ■ a ■ a/ (15) 

zrr sm d sm p 

S= max L ' J 

Proof. We start with the following addition theorem for the associated Legendre function of 



the second kind (see [3]), namely for 9 G (0,7r), 



00 

X 

n=m, 



where x > 1 is given by (15). By (8) and (11) the infinite series is uniformly convergent for 
9 G [0, vr]. Therefore, if we multiply both sides of (16) by sin0P™(cos0), where n' G No and 
integrate over 9 G (0, vr) we obtain (14). ■ 
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Corollary 1. Let n,m gNq with n > m, r,r' £ (0,oo), r / r' , 9' € (0,7r). Then 



/•TT 

Jo 



Q„^-l/2(x)P;^(cos0)^4i^d0 = ?|^^P-(cos0') 



where x > ^ is given by (15). 

Proof. Substitute i/ = -1 in (14) and use [12, (14.5.17)]. ■ 
Theorem 6. Let m,n G No with n>m, a,a' £ (0,oo), 9' G (0,-7r). Then 

Qm-i/2ix)Kicos9)V^d9 = 2^(-l)™^^^^^ 
X Vsinh a sinh a' sin 9' P;^ (cos 0') (cosh a< ) (cosh cj> ) , (17) 

where 

cosh^ a + cosh^ a' — sin"^ 9 — sin^ 9' — 2 cosh a cosh a' cos 9 cos 9' 
2 sinh C7 sinh a' sin sin 0' 

Proof. We start with the following addition theorem for the associated Legendre function of 
the second kind (see [6, (37)]), namely 



Qm-i/2{x) = vr(-l)'"\/sinh a sinh a' sin 9 sin 9' ^ (2n + 1) 



(n — m)!"' ^ 



(n + m)\ 

X P;^(cos0)P™(cose')^r(coshCT<)Q™(coshfT>), (19) 

where x ^ 1 is defined by (18). Note that x = 1 only if cr = cr' and 9 = 9', and in that case 
Qm-i/2ix) has a logarithmic singularity. If it 7^ cr' then (8), (10), (11) show that the series 
in (19) is uniformly convergent for 9 £ [0,7r]. Therefore, if we multiply both sides of (19) by 
Vsin P^ (cos 0) and integrate over 9 G [0,7r], then by (7) we have obtained (17). If a = a' then 
one may use (9), (10), (11) and the orthogonality relation (7) to show that the series in (19) (as 
a series of functions in the variable 9 S (0,7r)) converges in L^(0,7r). That is 

E ; ; '^»(^O^r(cosha<)Q^(cosha>) < 00, (20) 

■"^ I in + m ! 

n=m ^ ^ ' ■' 

where (pn £ -^^(0,vr) forms an orthonormal basis and is defined as 



for n = m,m + 1, . . . . Then by the asymptotics of and Q"^ (cf. (10), (11)) 

(cosh a<)g™ (cosh (j>) = O (n^™-^) as n ^ 00. 
Also from the estimate of P™ (9), (j)n{9') = O (1) . Therefore, 
(n — m)! 



(n + m)! 



6„(e')^"(cosha<)Q™(cosha>) = 0{n-^), 



and this implies (20) because ^ ^ < 00. Therefore, we again obtain (17). 

n=i 
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Theorem 7. Let m, n G No with < m < n, a,a' e (0, oo), 6*' G [0, vr]. Then 

r Qrn-i/2{x)Kicose)V^de = 2^i(-i)-^I^^^ 

Jo (n + mj! 

X Vcosh a cosh a' sin (cos ^O-CC^ sinh a<)Q™(i sinh a>) , 



where 



X 



sinh^ cj + sinh^ d' + sin^ 9 + sin^ 9' — 2 sinh cr sinh o"' cos cos 9' 



(21) 



(22) 



2 cosh a cosh u' sin 9 sin 0' 
Proof. We start with oblate spheroidal coordinates on M'^, namely 

X = a cosh a sin cos i;^, y = a cosh a sin sin (f), z = a sinh cr cos 9, 

where a > 0, a £ [0,oo), 9 £ [0,7r], cj) E [0,2-k). The reciprocal distance between two points 
X, x' € expanded in terms of the separable harmonics in this coordinate system is given in 
[9, (41), p. 218], namely 



1 



X — X 



£(2n + l)5^(-l)-6. 



n=0 



m=0 



(n — m)\ 
(n + m)\ 



cos{m{(j) — (j)')) 



X (cos 0)P;^ (cos 9')P^{i sinh cj< )Q;^ (i sinh cj> ) , 

where Cm '■= 2 — 5m,o is the Neumann factor [11, p. 744] commonly occurring in Fourier cosine 
series, with a' G [0, oo), 9' £ [0,7r], (j)' G [0,27r). Note that the corresponding expression given 
in [6, § 5.2] is given incorrectly (see [4]). By reversing the order of summations in the above 
expression and comparing with the Fourier cosine expansion for the reciprocal distance between 
two points, namely 



1 



1 



X — X' 



Tra 



V cosh a cosh a' sin 9 sin 9 



= ^ emcos(m(0 - (l)'))Qm^i/2{x), 



m=0 



where x > 1 is given by (22), we obtain the following addition theorem for the associated 
Legendre function of the second kind 



Qm-i/2{x) = ivr(— 1)"^^ cosh a cosh a' sin 9 sin 9' (2?! + 1) 



(n — m)\ 



(n + m)\ 



X P;^ (cos 0)P™ (cos e')^r {i sinh a< ) Q™ {i sinh a> ) . 



(23) 



If we multiply both sides of (23) by \/sin0P^(cos 9) and integrate over 9 G [0, vr], then by (7) we 
have obtained (21). We justify the interchange of integral and infinite sum as before by using 
the asymptotic formulas (12), (13). ■ 



Theorem 8. Let m, n G No with < m < n, a,a' £ (0, oo), 9' £ (0, vr). Then 
Q^^y2{x)Kicos9)V^9d9 = '^-^^lF^icos9')e-(-^'/'^(^>-'^<\ 
where if we define s = cosho", s' = cosho"', r = cos 9, r' = cos 9', then 



(24) 







X 



sm" »{g — t'Y + sin^ 9' {s — r)^ + [(s' — r') sinhcj — (s — r) sinhcr'] 



2 sin 6* sin 6" (s-r)(s'-r') 



(25) 
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Proof. We start with bispherical coordinates on M'^, namely 

a sin 6 cos 6 a sin 6 sin 6 a sinh a 



cosh a — cos 9 cosh a — cos cosh a — cos ^ 

where a > 0, a £ [0,oo), 6 G [0,7r], (p £ [0,27r). The reciprocal distance between two points 
X, x' G M'^ expanded in terms of the separable harmonics in this coordinate system is given in 
[9, (9), p. 222], namely 



1 °° 
= - V(cosha - cos 0) (cosh cj' - cos X] ^'^"^^^^^^"^"'"^^ 

(n — m)! 



a 

n=0 



^E^"^ ^ : p;r(cosg)p;r(cosgOcos(m(</>-<AO), 

where a' G [0,oo), 9' G [0,7r], i?!)' G [0,27r). By reversing the order of summations in the above 
expression and comparing with the Fourier cosine expansion for the reciprocal distance between 
two points, namely 

1 sJ (cosh a — cos 9) (cosh a' — cos 9') 

T| 7TT = / ■ ^ . ^, em COs(m(0 - ))Qrn-l/2{X), 

x-x' vra V sm 6* sm 6*' ^„ 

where x ^ 1 is given by (25), we obtain the following addition theorem for the associated 
Legendre function of the second kind 

Qm-1/2{X) = ^Vsin^sin^' p-^e-(-+i/2)(->-<)p-(cos0)P™(cose')- (26) 
' ^-^ (n + ml! 

n=m ^ ' 

\i a = g' and 9 = 9', then x = 1; and Qm-i/2ix) has a logarithmic singularity. Note that the 
corresponding expression given in [6, § 6.1, (45)] is given incorrectly (see [4]). If we multiply both 
sides of (26) by \/sin0P™(cos^) and integrate over 9 G [0,7r], then by (7) we have obtained (24). 
We justify the interchange of integral and infinite sum in the same way as in the proof of 
Theorem 6. ■ 

3.2 Order orthogonality for associated Legendre functions 
with integer degree and order 

In this subsection we take advantage of the order orthogonality relation for the Ferrers function 
of the first kind with integer degree and order (cf. [12, (14.17.8)]) 

Y>^{zos9)V^ {cos9)—d9 = - ^ (27) 
Q sm9 m [n — m)\ 

with m > 1. 

Theorem 9. Let m eN, n e No with 1 < m < n, 9' e [0,tt], (p, cf)' G [0, 27r). Then 

Pnicos-f)F'r {cos 9) ^d9 = — P™ (cos e')cos(m((^ -(/)')), 
Q sin 9 m 

where 

cos 7 = cos 9 cos 9' + sin 9 sin 9' cos{(j) — </>')■ 
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Proof. We start with the addition theorem for spherical harmonics (cf. [12, (14.18.1)]), namely 

P„(cos7)= j^-^,K{cose)F^{cos9')e^"^('f>-'f''\ (28) 

^-^ (n + ml! 

m=—n ^ ' 



where P„ : C — )• C, for n E No, is the Legendre polynomial which can be defined in terms of the 
terminating Gauss hypergeometric series (see for instance [12, Chapters 15, 18]) as follows 

— n, n + 1 1 — 

1 '^^y ■ 

We then take advantage of the order orthogonality relation for the Ferrers functions of the first 
kind with integer degree and order. If we multiply both sides of (28) by (sin0)~^ P™ (cos0) and 
integrate over Q G (0,7r), by using (27) we obtain the desired result. ■ 

Theorem 9, originating from (28), is the only example of a definite integral that we could find 
using the order orthogonality relation for the Ferrers functions of the first kind (27). Therefore 
we highly suspect that this result is previously known, and include it mainly for completeness 
sake. It would however be very interesting to find another example using this orthogonality 
relation. 



iFi 



3.3 Orthogonality for Chebyshev polynomials of the first kind 

Here we take advantage of orthogonality from Chebyshev polynomials of the first kind (cf. [12, 
§ 18.3]) 

TT 

Tm{cos9)Tn{cos6)d0 = —6m,n, (29) 

where : C — )• C, for n G No, is the Chebyshev polynomial of the first kind which can be 
defined in terms of the terminating Gauss hypergeometric series (see [12, Chapter 18]) 

Tn[z) = 2Fi{ ^ ; I . 

The Chebyshev polynomials of the first kind satisfy the identity [12, (18.5.1)] 

Tn{cos6) = cos{n6). 
Theorem 10. Let m,n £Z, a, a' £ (O,cxo). Then 



/ Qm-i/2{x) cos{n'>p)dtp = 7r(— 1)'"\/ sinh a sinh a' 
Jo 



2 

where 



F (tT/ — TTl -|- — 1 

p ^ ^1 Pn- 1/2 (cosh a< ) (cosh CT> ) , (30) 

1 (n + m + - ' ' ' 



X = coth a coth a' — csch a csch a' cos ^. (31) 

Proof. We start with toroidal coordinates on M'^, namely 

a sinh a cos (p a sinh a sin (p a sin tp 



cosh a — cos ip ' cosh a — cos ip ' cosh a — cos ip ' 
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where a > 0, a £ (0,oo), ip,4> £ [0,27r). The reciprocal distance between two points x, x' G 
is given algebraically by 



1 1 /(cosher — cos ^/^)( cosher' — cos V'') 



a V 2 sinh a sinh a' 

cosh a cosh a' — cos{^ — ip') 
sinh a sinh a' 



cos 



-1/2 



where {a' ,1^' , (p') are the toroidal coordinates corresponding to the point x'. Using Heine's 
reciprocal square root identity (see for instance [5, (3.11)]) 

^m.Qm-l/2{z)Tm{x) , 



m=0 

where z > 1 and x S [—1,1], we can obtain a Fourier cosine series representation for the 
reciprocal distance between two points in toroidal coordinates on M^, namely 



1 1 / (cosh fj — cos -i/;) (cosh fj' — cost/'') i/w^ 



llx — x'll vra V sinh £7 sinh cr' 

" " m=0 

where x > 1 is given by (31). We can further expand the associated Legendre function of the 
second kind using the following addition theorem (cf. [7, (8.795.2)]) 

oo 

Qm-i/2(x) = (-1)"'V sinh (7 sinh a' £„ cos{n{ip - ip')) 

n=0 

^ — T(^r-i/2(coshcT<)Q- i/2(cosha>). (32) 
r (n + m + |j ' ' 

Note that with the above addition theorem, we have the expansion of the reciprocal distance 
between two points in terms of the separable harmonics in toroidal coordinates 



1 



^ oo 

= — \/ (cosh a — cos ■0) (cosh a' — cos tp') (— l)™em, cos(m((/) — (j)')) 
llx — X' II vra 

" " m=0 

X cos(n(^ - ^')) ^ r ~ ^ ? ^r-i/2(cosh a<)g-_,/,(cosh a>) 

n=o i (^n + m + 2j 

(see also [6, § 6.2] and [4]). If we relabel ip — ip' ip and multiply both sides of (32) by cos{nip) 
and integrate over ^p G [0,tt], then by (29) we have obtained (30). The interchange of infinite 
sum and integral is justified by (10), (11). ■ 
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